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ABSTRACT 


Seismic  source  identification  work  during  this  reporting  period 
has  involved  the  development  of  several  new  analysis  tech¬ 
niques  for  almost  completely  automatic  production  measure¬ 
ment  on  seismograms,  and  for  less  routine,  more  detailed 
studies.  Stud'es  of  array  signal  detectability  and  sig.ial  and 
!  Oise  spatial  character  continue  using  data  from  both  Montana 
(LASA)  and  Norway  (NORSAR). 
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SUMMARY 


This  is  the  ninth  Semiannual  Technical  Summary  of  Lincoln  Laboratory's  work  for  the 
Advanced  Research  Projects  Agency  on  the  seismic  discrimination  problem  (Vela  Uniform). 

Identification  studies  using  short-period  and  long-period  small  and  large  arrays  are  con¬ 
tinuing.  Means  have  been  developed  for  applying  the  available  simple  waveform  discriminant 
parameters  to  large  numbers  of  events  with  a  high  degree  of  automaticity  (Sec.  1-A);  to  aid  the 
search  for  new  discriminants,  a  set  of  signal  processing  tools  which  provide  a  detailed  dissec¬ 
tion  of  individual  events  has  been  developed  (Sec.  I-B). 

Work  on  detection  and  location  at  single  stations  has  proceeded  in  several  directions.  LASA 
detection  and  two  methods  of  LASA  location  were  studied  using  data  recorded  during  the  Kuriles 
ocean  bottom  seismic  experiment  of  late  1966  (Sec.  11-A).  A  new  processing  method  for  high- 
resolution  space-frequency  analysis  of  signals  or  noise  (Sec.  II-B)  has  been  developed;  using 
this  and  other  tools,  a  study  of  A  .ntana  long-period  noise  was  completed  (Sec.  II-C).  Recent 
studies  of  signal  and  noise  properties  of  the  Norway  large  array  (NORSAR)  site  are  reported  in 
Sec.  li-D.  Section  II-E  summarizes  the  last  three  years'  work  on  measuring  and  using  LASA 
delay  station  corrections.  A  survey  of  work  elsewhere  on  signal-to-noise  improvement  by  seis¬ 
mometer  emplacement  in  deep  boreholes  has  been  made  (Sec.  II-F). 

A' description  of  the  latest  work  on  the  Seismic  Data  Analysis  Console  concludes  this  report. 

P.  E.  Green 
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SEISMIC  DISCRIMINATION 


I.  IDENTIFICATION 

A.  ARRAY  IDENTIFICATION  STUDIES  USING  LARGE  EVENT  POPULATIONS 

i  2 

Reports  published  earlier  have  described  the  application  of  long-period''  and  short-period 
discriminants  to  LASA  recordings  of  numbers  of  blast  and  earthquake  events.  The  data  base 
for  those  studies  was  oLiained  using  two  different  computers,  several  different  programs  and 
multiple  computer  runs  for  each  event  In  order  to  continue  such  large  population  experiments 
and  thus  gain  increased  understanding  of  the  existing  discriminants,  it  was  decided  that  the  pro¬ 
cedure  should  be  significantly  modified  and  simplified.  Therefore,  a  single  program  (Short- 
Period  Long-Period,  or  SPLP)  has  been  written  which  will  accomplish  allroutine  long-  and  short- 
period  analyses  of  array  station  (e.g.,  LASA)  data  in  a  single  pass.  The  SPLP  program  also 
stacks  data  on  a  composite  tape  for  future  reference. 

Figure  1  is  a  flow  chart  indicating  all  of  the  major  operations  done  by  the  program.  Except 
for  the  entry  of  seismograms  into  core,  no  input  or  output  has  been  shown.  With  respect  to  out¬ 
put  it  is  sufficient  to  say  that  all  measured  or  coinputeo  event  parameters  as  well  as  beam  wave¬ 
forms  and  spectral  curves  are  output  graphically  using  an  SC  4020  hard  copy  display,  and  are 
also  stored  on  a  composite  tape.  Although  many  options  are  available,  only  a  small  number  of 
input  parameters  is  changed  for  the  routine  processing  of  events.  In  particular,  the  program 
requires  rough  ..stimates  of  P  arrf  'al  time,  azimuih,  and  horizontal  phase  velocity  for  the  event. 
Arrival  times  within  20  seconds  are  quite  sufficient,  and  the  estimat.d  azimuth  (/3)  and  horizontal 
(v)  phase  velocity  supplied  to  the  program  can  be  quite  crude.  When  the  quality  of  0  and  v  are 
thought  to  be  particularly  poor,  a  larger  grid  of  SP  beams  can  be  requested.  Included  in  the  out¬ 
puts  stored  on  tape  are  the  unsmoothed  amplitude  spectral  components  of  the  first  20  seconds  of 
the  P-phase  which  arc  available  for  use  in  deriving  possible  improved  versions  of  the  spectral 
ratio  discriminant. 

The  routine  processing  by  SPLP  differs  in  several  ways,  in  addition  to  the  obvious  fact  that 
much  is  automatic,  from  the  processing  done  for  earlier  discrimination  studies.  The  three  most 
significant  differences  are  the  following.  First,  bandpass  filtering  is  done  by  Fourie**  transform¬ 
ing,  shaping,  and  inverse  transforming  rather  than  recursive  or  time  domain  convolutional  fil¬ 
tering.  Second,  SP  spectra  are  obtained  by  Fourier  transforming  and  spectral  ratio  is  obtained 
from  the  transformed  data.  Lastly,  chirp  filtering  is  done  by  applying  a  quadratic  phase  shift  to 
Fourier  transformed  data  and  then  inverse  transforming  rather  than  by  time  domain  convolution. 

A  population  of  events  detected  and  recorded  at  LASA  has  been  compiled  for  analysis  using 


*  SPLP  will  accept  either  Lo-Rate  or  Hi-Rate  tape  from  LASA  as  input.  If  Lo-Rate  tape  is 
supplied,  then  both  long-period  and  short-period  analysis  is  done.  Otherwise,  only  the  short- 
period  analysis  is  completed.  The  program  requires  about  20  minutes  of  IBM  360/67  time  to 
completely  process  an  event.  This  is  reduced  to  15  minutes  if  only  short-period  data  are  avail¬ 
able.  The  program  has  been  made  sufficiently  general  so  that  it  can  be  applied  to  short-period 
data  recorded  at  other  array  sites  with  a  minimum  of  effort. 
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Section  I 


ENTER  UF  TO  102  SECONDS  OF  SF  DATA 
FROM  UF  TO  21  CHANNEIS  INTO  COM 
MEMOEY 


CHECK  DATA  OUAUTY  AND 
HJECT  lAO  DATA 


lANDFASS  FILn*  DATA 

0.6  -  2.0  Ht 


FORM  lEAMS  AT  AND  AROUND 
NOMINAL  V/AVENUMUR  OF  EVENT 
OYFKALLY  100  -  400  REAMS) 

STATION  CORMCTIONS  USED 


SKMINE  INIERVAL  IN  WHICH  I 
SEARCH  FOR  LF  RAYIEIOK  WAVES! 


MFEAT  FOR 
EACH  OF  3 
lONG-FERlOO  < 
COMFONENTS 
(Z,  N-S,  E-W) 


ENTER  UF  TO  34  MINUTES  OF  ONE 
COMFONENT  OF  LF  DATA  FROM 
UF  TO  21  SITES 


CHECK  DATA  QUAUTY  AND 
REJECT  UD  DATA 


RANDFASS  FILTER  DATA 
0.02S -O.OSS  Ha 


SAVE  WAVENUMRER  OF  REAM 
WITH  URGEST  FEAK-TO-FEAK 
OUTFUT 


FORM  REAM  STEEMD 
TO  NOMINAL  RAYLEIOH 
VELOCITY  OF  3. 7  km/MC 


USING  REST  FILTEMD  REAM  AND 
ITS  WAVENUMRER  COMFUTEi 
FERIOD 
AMFUTUDE 
VELOCITY 
AZIMUTH 
DISTANCE 
UTITUDE 
lONonUOE 

ROOY-WAVE  MAGNITUDE 
COMFUlXITY 

AND  CHECK  !>0(  FOSSIRIE  DEPTH 
I'HASES. 


COMPUTE  HORIZONTAL 
SEISMOGRAMS  RADIAL 
AND  TRANSVERSE  TO 
DIRECTION  OF  ARRIVAL 
OF  EVENT 


APFIY  FAMILY  OF  CHIRP 
FILTERS  TO  LPZ  REAM 
AND  SAVE  THAT  WITH 
largest  ZERO-TO-FEAK 
VALUE 


FORM  UNFILHMD  ream 
CORRESPOF4Dlt4G  TO  REST 
FILTERED  MAM 


COMFUTE  SURFAQ  WAVE  MAGNITUDE 
FROM  CHIEF  FILnUED  LPZ  REAM 


COMFUTE  SFECTRUM  OF 
UNFILTEMD  REAM  AND 
COMFUTE  SFEaRAl  RATIO 


APPLY  CHIRP  FILTER  TO 
RADIAL  REAM 


Fig.  1.  Sequence  of  outamatic  seismogrom  reading  operations  in  SPLP  program. 
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Section  I 


SPLP.  It  includes  some  30  tclescismic  explosions,  35  large  shallow  Sino-Soviet  events  before 
1968,  and  LASA-detected  ;.mo-Soviet  events  between  i  January  1968  and  8  May  1968.  So  far 
SPLP  has  been  successfully  applied  to  more  than  120  of  the  175  events  of  this  population  tiiat  we 
intend  to  process.  R.  T.  Lacoss 

B.  SPECIAL  METHODS  FOR  DETAILED  ANALYSIS  OF  INDIVIDUAL  EVENTS 

The  procedures  described  in  Sec.  1-A  are  being  used  to  measure  a  few  simple  parameters 
on  a  large  number  of  events  for  statistical  evaluation.  The  "console"  development,  described 
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earlier  and  in  Sec.  IV-A,  is  motivated  by  a  similar  philosophy,  namely  the  carrying  out  of  a 
number  of  pre-programmed  measurements  on  recorded  data.  In  this  Section  we  describe  a  se¬ 
ries  of  programs  based  upon  the  contrasting  app-oach  of  performing  several  very  detailed  an¬ 
alyses  on  the  LASA  waveforms  of  a  small  number  oi  events  in  tlie  nope  of  discovering  details  of 
the  wave  motion,  as  seen  by  the  array,  which  may  be  made  the  basis  of  new  discriminants.  The 
emphasis  is  on  providing  useful  hard  copy  displays,  in  addition  to  conventional  seismograms, 
to  be  studied  by  a  human  observer.  Three  such  displays,  which  are  now  in  operation,  will  be 
described.  Tliey  are  all  new  and  improved  versions  of  techniques  which  have  been  described  in 
previous  reports. 

The  first  is  a  velocity  spectral  analysis,  called  VESPA,  which  accepts  21  subarray  wave¬ 
forms  and  outputs  signal  power  as  a  function  of  wave  velocity  and  time.  In  this  program,  which 
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replaces  the  "beam  power  display"  described  earlier,  '61  beams  are  formed  from  the  bandpass 
filtered  (0.6  to  2.0Hz)  subarray  outputs,  which  may  be  weighted  for  control  of  array  response. 

The  beams  are  aimed  at  the  bearing  of  the  event  bcbig  studied,  at  pre-assigned  horizontal  phase 
velocities  corresponding  to  uniform  increments  in  slowness  of  0.2  second  per  degree,  from  0  to 
12  seconds  per  degree.  Moveouts  are  computed  from  the  1968  Herr  in-Taggart  travel  times  and 
the  best  available  station  corrections  are  used.  The  beam  output.?  are  formed  once  each  LASA 
time  sample  (i.e.,  at  20  Hz)  and  its  squares  are  summed  over  a  1 -second  time  irtt-rval  (integrate 
and  dump)  to  produce  output  powers  at  a  1-Hz  rate.  These  powers  arc  summed  between  pre- 
assigned  limits  (say  20  seconds  covering  the  main  P-arrival)  and  displayed  as  a  spectral  density 
curve,  i.e.,  power  as  a  function  of  slowness.  In  addition,  the  61  beam  powers  are  treated  as 
functions  of  time  over  one  or  more  1-minute  framer  and  displayed  as  contour  diagrams  of  power 
vs  slownes'i  and  time.  The  contours  are  curves  oi  constant  numbers  of  cib  down  from  the  peak 
value  observed  in  the  first  minute,  and  the  display  show.s  the  detailed  variaiion  with  time  of  the 
distribution  of  signal  power  with  re.'^pect  to  slowness. 

Ir  Fig.  2(a)  we  see  such  a  2-minute  frame  for  an  East  Kazakhstan  presumed  underground  test. 
The  LASA  D,  K,  and  F  rings  were  used,  with  weights  applied  to  equalize  the  peak-to-peak  ampli¬ 
tudes.  Three  of  the  constituent  filter  jd  seismograms  are  shown  in  Fig.  2(h).  The  total  spectrum 
of  power  vs  slowness  for  the  20-second  interval  marked  on  Fig.  2(a)  is  given  in  Fig.  2(c).  The 
main  lobe  agrees  very  well,  down  to  about  lOd'-,  ivith  a  theoretical  pattern  for  a  wide-band 
source,  monochromatic  in  velocity.  Analogous  display.?  are  shown  in  h'igs.  3(a),  (b),  and  (c)  for 
a  presumed  underground  test  at  Novayn  licmlya,  using  the  same  collection  of  12  I.A.SA  sen.sors. 
The  event  is  much  more  complex,  and  signal  energy  with  very  nearly  the  I'-vclocity  can  be  seen 
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Section  I 


to  arrive  in  substantial  amounts  for  over  1  minute,  while  the  first  20  seuon.ls  still  show  a  rea¬ 
sonable  velocity  spectrum.  The  phases  PcP  and  PP  are  also  visible  in  the  temporal  display, 
which  in  this  case  shows  3  minutes  of  data. 

The  second  display  is  a  new  version  of  the  previously  described  "sonogram"^  program.  A 
single  input  waveform  is  pas.sed  through  the  digital  equivalent  of  a  bank  of  31  narrow-band  (three- 
pole  Dutterworth)  filters,  all  having  the  same  Q.  The  power  output  of  each  filter  is  summed  for 
1  second  (20  samples)  and  these  data  are  displayed  as  contour  diagrams  of  power  vs  frequency 
and  time  for  1 -minute  frames,  exactly  as  in  the  VESPA  program.  The  program  actually  gener¬ 
ates  sonograms  for  up  to  live  input  channels  .simultaneously,  and  as  usually  applied,  .one  channel 
is  a  preformed  beam  steered  at  the  event.  An  example  of  the  sonogram  display,  using  the  same 
Novayr  'lemlya  event  of  Fig.  3,  is  shown  in  Fig.  4.  The  characteristic  delay  in  the  arrival  of 
lower  frequency  energy  is  due  to  the  fact  that  all  filters  have  the  same  Q  (bandwidth  is  10  percent 
of  center  frequency),  hence  the  low-frequency  filters  have  greater  delay  time.  It  has  been  ob¬ 
served  that  the  frequency  of  peak  energy  varies  considerably  from  subarray  to  subarray  for  a 
single  event,  and  is  only  weakly  correlated  with  the  dominant  period,  as  read  from  the  original 
seismograms. 

The  third  display  is  a  particle  motion  analysis,^  using  the  short-period  triaxial  sensors  in¬ 
stalled  in  subarray  '‘2.  The  horizontal  records  are  conve*  ted  to  radial  and  transverse  compo¬ 
nents  and  particle  motion  tracings  of  vertical  .'s  radial  and  radial  vs  transverse  are  generated 
for  a  sequence  of  2-second  frames.  The  expected  angle  of  particle  motion  in  the  vertical-radial 
plane  is  displayed,  along  with  raw  reismograms  and  a  time-smoother  version  of  the  squared  ver¬ 
tical  trace.  This  last  waveform  is  roughly  similar  to  the  AWRE  correlogram^  output. 

None  of  these  programs  has  as  yet  been  run  on  significant  numbers  of  events,  but  it  is  our 
intent  to  use  th^m  for  detailed  study  of  representative  earthquakes  and  explosions  in  the  hope  of 
finding  new  P-wavoform  discriminants  and  of  elucidating  the  structure  of  the  P-coda. 

E.  J.  Kelly 
L.  T.  Fleck 
P.  E.  Green 
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n.  SINGLE-STATION  DETECTION  AND  LOCATION 


A.  COMPARISON  OF  LASA  OBSERVATIONS  WITH  KURILES  OCEAN  BOTTOM 
NET  OBSERVATORIES 

The  Kurile  Islands  offshore  experiment  condu-jted  at  the  end  of  1966  provided  the  opportunity 
of  obtaining  accurate  locations  of  several  small  underwater  explosions  and  many  small  local  and 
near  regional  earthquakes  from  a  seismically  active  area.  This  population  of  events  has  been 
used  to  better  determine  the  Montana  LASA  beam  detection  threshold  for  teleseismic  P  waves 
and  the  epicenter  location  error  by  the  use  of  both  plane  wave  approximation  to  the  arrival  times 
at  selected  subarrays  and  beamsplitting.  A  detailed  report  on  this  work  is  in  preparation. 

The  beam  detection  threshold  for  the  Kuriles  was  obtained  from  42  events  thus  far  processed. 
If  a  seismicity  curve  with  a  unity  slope  is  used,  the  75  percent  curuulative  detection  threshold  is 
3.8.  If  a  curve  is  fit  to  the  data,  a  slope  of  0.75  is  obtained  with  a  75  percent  cutnulative  detec¬ 
tion  threshold  of  3.65.  Figure  5  summarizes  the  results. 

The  accuracy  of  epicenter  location  by  either  plane  wave  approximation  (fitting  a  plane  wave 
to  individual  subarray  times)  or  beamsplitting  is  strongly  affected  by  station  time  errors.  Station 

g 

corrections  obtained  from  Teledyne  were  used  to  compensate  for  these  time  anomalies.  The 
plane  wave  epicenter  locations  produced  epicenter  errors  shown  in  Fig.  6(a),  while  the  epicenter 
location  errors  from  the  beamsplitting  method  are  shown  in  Fig.  6(b).  (See  also  Fig.  7  of  the 
31  December  1967  Seismic  Discrimination  Semiannual  Technical  Summary  Report.)  The  epi¬ 
centers  obtained  from  the  plane  wave  approximation  method  were  computed  to  0.1  degree;  how¬ 
ever.  the  epicenters  obtained  from  the  beamsplitting  method  are  quantized  to  0.5-degree  incre¬ 
ments  of  latitude  and  longitude  such  that  the  best  beam  must  fall  on  one  of  these  increments.  The 
epicemer  location  error  plotted  for  beamsplitting  can  be  considered  to  be  the  upper  bound  of  the 
location  error  because  of  the  coarseness  of  the  increments. 

The  data  base  used  to  get  the  above  detection  and  location  results  contains  84  events  that 

occurred  in  the  Kurile  Inlands  region  between  23  October  and  12  December  1966.  The  epicenter? 

were  obtained  from  four  sources:  (1)  the  Coast  and  Geodetic  Survey  PDE  list.  (2)  Texas  Instru- 

g 

ments,  Inc.,  analysis  of  ocean  bottom  seismometer  (OBS)  data,  (3)  Lincoln  Laboratory  analysis 
of  OBS  data,  and  (4)  the  LASA  station  bulletin.  The  LASA  station  bulletin  epicenters  were  con¬ 
sidered  the  least  accurate  not  only  because  of  the  limited  baseline  of  the  array,  but  also  because 
the  LASA  epicenter  location  computer  program  had  not  been  fully  implemented  with  station  cor¬ 
rections  at  that  time. 

A  further  discussion  of  category  (3)  is  in  order.  In  order  to  increase  the  number  of  small 
events  contributed  by  the  OBS  data,  epicenters  were  computed  for  sets  of  events  not  located  by 
Texas  Instruments.  In  order  tcj  locate  these  events,  we  assume  a  surface  focus  (Texas  Instru¬ 
ments'  location  included  a  depth  determination  and.  therefore,  required  detection  at  more  of 
the  OBS  instruments)  so  that  only  two  P  arrivals  and  ot.e  S  arrival  were  required  to  determine 
an  epicenter.**^  If  only  two  P  arrival?,  and  one  S  arrival  are  used,  two  epicenters  are  produced. 
Any  additional  readings  serve  to  resolve  the  ambiguity  and  refine  the  epicenter.  The  accuracy 
of  the  resulting  epicenters  depends  ujjon  good  P  .md  S  readings  and  upon  a  knowledge  of  the  P 
speed  and  the  S-P  interval  as  a  function  of  distui.cc.  We  first  attempted  to  get  S-P  intervals 
using  the  underwater  explosions  for  which  an  accurate  epicenter  was  known,  but  the  resulting 
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LASA  BODY  WAVE  MAON'TUdE 


(a)  LASA  epicenters  derived  from  plane  wave  fit  to  individual 
subarray  times  ("triangulatian"). 
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LAM  BOOY  WAVE  MASNITUOE 

(b)  LASA  epicenters  derived  from  beamsplitting  (determining 
beam  with  largest  autput). 


Fig.  6.  Disagreement  between  LASA  epicenters  and  epicenters  from  either  USCGS 
ar  Kuriles  acean  battam  seismometer  net. 
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S-P  intervals  showed  a  large  scatter,  probably  because  of  the  water  part  of  the  propagation  path. 
Tne  population  of  earthquakes  located  by  the  USCCS  produced  more  consistent  S-P  intervals  and 
enabled  a  good  determination  of  the  S-P  separation  with  distance.  Figure  7  shows  the  data  and 
Ue  resulting  curve.  ^  Sheppard 

B.  HIGH-RESOLUTIOK  FREQUENCY-WAVENUMLER  SPECTRUM  ANALYSIS 

The  frequency-wavenumber  spectrum  of  seismic  noise  orovides  the  information  about  the 
distribution  of  noise  power  with  frequency  and  fo.  a  fixed  frequency  reveals  the  velocity  and 
direction  of  propagation  of  the  noise.  A  high-resolution  method  for  measuring  the  frequency- 
wavenumber  spectrum  of  seismic  noise  has  been  developed.  A  detailed  description  of  this  pro¬ 
gram  as  well  as  the  conventional  method,  have  been  given  in  a  recent  report.^  ^ 

The  conventional  method  of  frequency-wavenumber  spectrum  analysis  is  based  on  the 
estimate 
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where  f  is  the  i^equency,  k  is  the  vector  wavenumber,  K  is  'he  number  of  sensors,  Xj  is  the 
vector  position  of  the  sensor  and  {r^(f)}  is  the  estimate  of  the  spectral  matrix  of  the  noise. 
The  high-resolution  estimate  is  defined  as 
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where  {qj^lf))  is  the  inverse  of  the  spectral  matrix  {f^^(f)).  It  can  be  shown  that  P'(f,  k)  is  the 
output  of  a  maximum-likelihood  filter,  whose  design  is  determined  by  the  noise  data  and  is 
different  for  each  wavenumber  k,  which  passes  undistorted  any  monochromatic  plane  wave 
travelling  at  a  velocity  corresponding  to  the  wavenumber  k,  and  suppresses  in  an  optimum 
manner  the  power  of  those  noise  waves  travelling  at  velocities  corresponding  to  wavenumbers 
other  than  k.  It  should  be  noted  that  the  amount  of  computation  required  to  obtain  P'  is 

A 

almost  the  same  as  that  to  get  P,  since  only  an  additional  Hermitian  matrix  inversion  is 
required. 

Two  examples  which  illustrate  the  results  obtained  with  both  the  conventional  and  high- 
rosolution  frequency-wavenumber  programs  are  given  in  Figs.  8(a)  and  (b).  In  both  rases  the 
entire  long-period  vertical  array  at  LASA  was  employed.  It  has  been  shown  that  the  noise  in 

the  0.025-  to  0.05-Hz  band  consists  of  a  component  propagating  as  fundamental-mode  Rayleigh 

il 

waves,  and  a  nonpropagating  or  incoherent  component.  Thus,  the  wavenumber  structure  of 
the  noise,  at  a  particular  frequency,  should  consist  of  an  arc,  or  arcs,  whose  extent  corre¬ 
sponds  to  the  range  of  the  azimuths  of  the  sources  of  the  propagating  component  c.  the  noise. 

In  Fig.  8(a)  the  noise  is  seen  to  consist  primarily  of  a  component  propagating  from  the 
North.  Hence,  as  expected,  the  conventional  wavenumber  results  appear  to  be  very  similar  to 
a  l.ASA  bcamforming  array  response  pattern.  The  high -resolution  results  delineate  the  wave¬ 
number  structure  as  a  point  source  more  clearly  than  does  the  conventional  method.  This  figure 


15 


Section  II 


f  > 

WAVCNUMaCR  (C|icl<l/knl  W/kVENUMtEO  (crclu/hni 


FREOUENCT  •  0.01  HI 

T  *p->  *T  NOiSE  SAMPLE 
21^1u:00  TO  C0'.10:00 

(a)  Point  ‘irce. 


Fig.  8.  Comparison  of  conventional  (left)  and  high  resolution  (right)  LASA  long-period  orray 
wavenumber  plots  for  two  interesting  noise  samples. 
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(b)  Single  velocity,  oil  azimuths. 


Fig.  8*  Continued. 
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shows  that  there  is  an  improvement  in  wavenumber  resolution  of  the  high-resolution  method 
relative  to  the  conventional  method  of  about  a  factor  of  four.  The  lesults  of  Fig.  8(b)  are  also 
interesting,  as  they  show  a  360-degree  azimuthal  sp.’-ead  for  the  wavenumber  structure  with  a 
variable  energy  density  along  this  circle.  In  addition  to  the  noise  analysis,  an  effort  is  also  in 
progress  to  apply  the  high-resolution  method  to  the  short-period  P-wave  data  in  an  atten'p:  to 
improve  the  location  capability  of  '?rge  arrays. 

J.  Capon 

C.  NOISE  MEASUREMENTS  IN  MONTANA 

The  long-period  noise  in  the  20-  to  40-second  period  .  ange  limits  the  identification  level  at 
which  the  surface-wave,  body-wave  discriminant  can  be  applied  at  LASA.  Therefore,  an  inves¬ 
tigation  was  made  to  determine  the  sources  and  properties  of  this  noise.  Only  the  long-period 
vertical  array  at  LASA  was  considered. 

The  conventional  and  high-resoluti''n’''  frequency-wavenumber  spectra  of  the  noise  were 
measured,  as  well  as  coherence.  These  data  shew  that  the  noise  consists  of  two  components. 

One  component  propagates  across  the  array  as  'undamental-mode  Rayleigh  waves  and  is  thought 
to  be  cau'  ;d  by  the  action  of  surf  on  coastlines.  The  other  component  is  nonpropagating  and  has 
been  found  to  be  incoherent  over  distances  greater  than  7.5  km. 

It  is  important  to  establish  the  origin  of  the  nonpropagating  seismic  noise  since  if  means 

could  be  found  to  suppress  it,  thus  exposing  the  propagating  com  'onent  more  prominently, 

maximum-likelihood  processing  could  be  used  effectively  against  the  latter  to  improve  LP  de- 

12 

tection  thresholds.  Previous  results  due  to  Haubrich  and  MacKenzie  indicate  that  the  non¬ 
propagating  noise  might  be  due  to  the  elastic  loading  on  the  ground  by  the  earth's  atmosphere. 
Thus,  an  effort  was  made  to  determine  the  coherence  between  the  seismic  noise  and  the  itmos- 
pheric  fluctuations  as  recorded  on  microbarograph  sensors  at  LASA. 

Two  examples  in  which  coherence  was  measured  at  site  A|^  a.  e  given  in  Figs.  9(a)  and  (b), 
where  Fig.  9(b)  is  more  typical  of  the  behavior  of  the  coherency  when  any  coherence  is  measured 
at  all.  In  Fig.  9(a)  the  coherence  at  0.03  Hz  is  quite  high,  about  0.6,  and  drops  to  about  0.5  .\nd 
0,4  at  0.04  and  0.05  Hz,  respectively.  The  ratio  of  nonpropagating  noise  to  the  total  noise,  for 
this  noise  sample,  was  measured  as  0.9  at  0.03,  0.04  and  0.05  Hz.  The  measured  coherence  is 
not  as  high  as  this  at  these  frequencies,  but  is  the  largest  that  was  ever  measured  in  the  17  noise 
samples  examined.  The  results  of  Fig.  9(b)  are  more  typical  in  that  whenever  coherei.ee  is  meas¬ 
ured,  it  tends  to  be  high  at  0.02  Hz  and  then  drops  to  the  level  for  incoherent  noise  at  0.03  to 
0.05  Hz.  The  ratio  of  nonpropags  ting  to  total  noise  for  the  noise  sample  used  in  Fig.  9(b)  was 
measured  as  0.60,  0.45  and  0.35  at  0.03,  0.04,  and  0.05  Hz,  respectively.  The  coherence  never 
reaches  a  level  compatible  with  the  amount  of  nonpropagating  noise.  Thus,  it  may  be  concluded 
that  there  is  relatively  little  coherence  between  the  seismic  noise  and  the  microbarograph  signals 
in  the  20-  to  40-second  period  range. 

An  attempt  was  made  to  correlate  the  amount  of  nonpropagating  seismic  noise  power  and  the 
amount  of  microbarograph  noise  power  in  the  20-  to  40-second  period  range.  These  quantities 
were  measured  at  site  A(J  for  17  noise  samples  and  the  results  are  shown  in  Fig.  10.  This  figure 


*See  Sec.  II-B. 
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shows  that  there  is  a  definite  trend  for  the  power  level  of  the  nonpropagating  noise  to  increase 
when  the  power  level  on  the  microbarograph  increases.  Thus,  this  is  evidence  that  the  non¬ 
propagating  seismic  m  ise  is  caused  by  the  elastic  loading  on  the  ground  by  the  earth's  atmos¬ 
phere.  The  nonpropagating  seismic  noise  may  be  caused  by  the  cumulative  effects  of  many 
independent  atmospheric  pressure  fluctuation  cells  within  a  radius  of  a  few  km  of  the  seismom¬ 
eter  and  this  may  be  the  reason  for  the  low  coherence  between  this  noise  waveform  and  a  single 
microbarograph  record  obtained  at  the  sa.ne  location. 

A  detailed  description  of  these  results  has  been  included  in  a  recent  report describing 
the  technique  discussed  in  the  preceding  section. 

J.  Capon 

D.  NORWAY  SIGNAL  AND  NOISE  PROPERTIES 

Analysis  of  noise  and  signal  data  from  the  Norwegian  Seismic  A»  ray  (NOP.SAR)  site  has 
continued.  The  earlier  short-period  results  ^  liave  been  substantiated  and  extended.  Some 
properties  of  the  long-period  noise  and  signal  have  also  been  obtained. 

The  effect  of  seismometer  burial  and  local  weather  conditions  upon  noise  levels  has  been 
thoroughly  investigated  using  data  from  the  Kaldalen  site.  Our  conclusion  is  that  instruments 
which  are  60  meters  into  bedrock  will  seldom,  if  ever,  be  more  than  3  db  quieter  in  the  signal 
band  than  instruments  only  2  meters  into  bedrock.  Average  and  peak  winds  were  obtained  from 
data  recorded  by  the  baldalen  weather  station  from  28  December  1967  to  22  March  1968.  Simul¬ 
taneous  spectra  were  obtained  from  instruments  at  2  and  60  meters  at  the  Faldalen  site.  The 
wind  conditions  during  these  measurements  ranged  from  calm  to  periods  with  gusts  to  45  mph 
and  sustained  winds  to  24  mph.  These  tatter  conditions  are  among  the  most  severe  during  the 
observation  period.  Figure  11  shows  the  range  of  noise  attenuation  that  was  recorded  at  the 
two  60-meter  holes  in  the  Faldalen  site.  It  can  be  seen  that  almost  3  db  of  attenuation  was  ob¬ 
tained  at  least  above  2.0  Hz  at  one  of  the  60-meter  holes.  However,  below  1.5  Hz  not  even  that 
small  amount  of  noise  attenuation  was  observed.  In  fact,  below  1.5  Hz  the  greatest  attenuation 
at  60  meters  was  only  about  1.5  db. 

it  has  been  verified  that  a  system  with  about  4  km  minimum  spacing,  such  as  the  0yer  sub¬ 
array.  will  result  in  noise  rejection  in  the  normal  signal  band  above  0.6Hz.  This  was  done 
by  making  use  of  frequency  wavenumber  analysis  of  the  noise  as  well  as  by  actually  forming 
beams  from  the  0yer  site.  This  work  as  well  as  further  frequency  wavenumber  analysis  of 
Faldalen  and  Ringsaker  data  has  indicated  that  surface  modes  with  horizontal  phase  velocities 
in  the  3  to  4  km/sec  range  tend  to  make  up  the  dominant  component  of  organized  short -period 
noise  in  Norw'ay. 

Extreme  amplitude  variations  over  the  (?!)yer  subarray  (approximately  20  km  aperture)  are 
approximately  a  factor  of  2  which  is  to  he  compared  to  a  factor  of  3  for  7-km  subarrays  in 
Montana. 

Figure  12  shows  the  typical  range  of  power  spectral  densities  observed  in  Norway.  The 
curves  are  uncorrected  for  instrument  response.  The  general  Norway  noise  level  has  been 
verified  using  more  than  20  spectra  from  the  Faldalen  noise  survey  site  and  a  smaller  number 
from  the  0yer  site.  .Although  the  Norway  noise  level  is  higher  at  1  Hz  than  it  is  in  .Montana, 
there  is  a  more  rapid  drop  with  increasing  frequency.  Short-period  P-magnitudes  at  0yer  have 
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Fig.  13.  Spectral  differences  in  P-pha$e  at  Nonway  and  Moi.i"  \a  in  the  form  af  o  histogram 
of  dominant  periods. 
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been  compared  with  USCGS  magnitudes  fur  25  events  and  it  was  found  that  the  averagi-  difference 
'er  ~  ‘•'O-l®  and  —0.09  for  events  from  the  east  (20°  -  90° )  and  west  (250°  -  310° ) 

respectively,  and  +0.07  averaged  overall  azimuths, 

One  of  the  most  notable  differences  between  signals  seen  .at  L.^SA  and  NOHSAU  seems  to  be 
the  wider  spread  in  the  frequency  content  of  the  NORSAR  signals,  Figures  13(a)  and  (b)  show 
histograms  of  the  dominant  periods  of  a  set  of  events  seen  at  NORSAR  and  LASA,  respectively. 
The  dominant  periods  at  LA.SA  cluster  closely  around  0,9-second  periods  while  the  dominant 
periods  at  NORSAR  spread  from  0,5  to  1,5  seconds. 

As  another  demonstration  of  the  wider  variation  of  spectral  characteristics  in  signals  seen 
at  NORSAR  compared  to  those  at  LASA,  Fig.  14(a)  and  (b)  show  the  spectra  of  an  event  in  .lapan 
as  seen  at  LASA  and  NORSAR.  For  this  event,  the  NORSAR  signal  has  significantly  wider  band¬ 
width  than  the  LASA  signal.  In  contrast,  Figs.  15(a)  and  (b)  show  the  spectra  of  an  event  in  the 
mid-Atlantic  ridge,  also  approximately  equidistant  from  the  two  sites,  for  which  the  spectra 
are  very  similar.  It  appears  that  the  NORSAR  essentially  provides  a  wider  band  receiver  than 
LASA  so  that  signal  characteristics  seen  at  NORSAR  may  be  more  influenced  by  the  actual  source 
function  than  at  LASA. 

A  number  of  power  spectral  densities  have  been  estimated  using  long-period  data.  The 
purpose  was  to  obtain  information  concerning  the  LP  noise  level  and  the  relatii  e  amounts  of 
noise  on  vertical  and  horizontal  instruments. 

Figures  16(a)  and  (b)  show  some  of  the  spectra  obtained  from  three  of  a  total  of  five  ob¬ 
servation  periods.  These  spectra  arc  not  corrected  for  instrument  response  and  are  relative  to 
1  (mp)^/Hz  at  0.04  Hz.  The  instruments  are  of  the  same  type  as  those  used  in  .Montana.  Other 
measurements  "/“re  made  for  14  .lanuary  and  2  February.  Figure  16(a)  sliows  S|)e<-tra  for  all 
three  components  on  28  January  which  was  the  observation  time  wliicli  gave  the  highest  noise 
level.  Assuming  a  signal  band  from  40-  to  20-second  periods  we  note  that  the  peak  in-band 
noise  is  up  to  55  db  above  the  reference  level  on  the  vertical  instrument  and  abu-.t  50  db  on  the 
horizontal  instruments.  Out  of  band,  the  peaks  for  the  16-  and  8-socond  microseisnis  are  up 
to  almost  60  db  on  verticals,  slightly  less  on  horizontal  instruments.  Figure  16(b)  shows  power 
spectral  densities  for  28  January  as  well  as  for  periods  about  one  day  earlier  and  latei’,  The 
in-band  power  is  reduced  from  a  few  to  10  db  on  those  days.  The  powei-  at  8  and  16  seconds  is 
reduced  even  more.  During  this  entire  3-day  period  there  was  a  developed  low  pi’cssure  area 
just  off  the  northwest  coast  of  Norway,  which  is  thought  to  be  the  source  of  the  disturbance, 

Our  observations  of  spectra  on  14  January  and  2  February  gave  noise  peaks  only  up  to  about 
40  db  above  the  reference.  A  number  of  observations  of  LP  noise  in  Montana  indicate  that  the 

typical  in-bar.d  vertical  noise  is  about  40  db  and  that  it  can  be  several  dl)  higlier  at  the  8  and  16 
11 

second  peaks.  The  horizontal  noise  in  Montana  is  often  up  to  10  dlt  more  intense  than  the 
vertical  noise. 

It  appears  from  the  above  that  the  vertical  long-period  noise  component  in  Norway  can  be 
as  quiet  as  typical  Montana  conditions.  However,  microscismic  storms  may  be  mure  intense 
in  Norway  resulting  in  vertical  noise  levels  10  to  20  db  higher  than  tluise  typical  in  .Montana. 
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Fig.  14.  Comparative  spectra  af  some  event  at  LASA 
and  NORSAR:  earthqu^e  aff  East  toost  of  Honshu, 
11  January  1968,  1612  Z,  mb(LASA)  =  4.8.  (o)  LASA 
beam,  (b)  NORSAR  0y#r  subarroy  beam. 


Fig.  15.  Comparative  spectra  af  some  event  at  LASA 
and  NORSAR;  earthquake  an  Mid-Atlantic  Ridge, 
8  January  1968,  2022  Z,  tnj,(LASA)  =  5.3.  (a)  LASA 
beam,  (b)  NORSAR  0yer  subarray  beam. 
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(a)  Vertical  component  before,  during  and  after 
micraseismic  noise  peak. 


(b)  Comparison  of  vertical  and  horizontal  output 
at  micraseismic  storm  peak. 


Fig.  16.  Norway  (Foldalen)  LP  noise  spectra. 
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Of  course,  very  ituet  -e  storms  can  similarly  raise  noise  levels  in  Montana  but  the  intensity  and 
duration  is  likely  to  be  greater  in  Norway.  .Apparently  the  horizontal  instruments  in  Norway 
have  .somewhat  less  noise  than  the  verticals.  This  is  particularly  true  at  periods  greater  than 
25  or  30  seconds  and  when  the  vertical  component  noise  level  is  high.  The  in-band  horizontal 
noise  level  in  Norway  is  typically  a  few  db  less  than  that  in  Montana  and  should  never  be  any 
worse  than  the  typical  Montana  level  even  during  intense  microseismic  storms  in  Norway. 

From  a  study  of  surface  waves  on  nine  events  seen  at  both  NORSAH  and  LASA,  it  was  found 
that  on  the  average,  values  at  the  two  stations  are  equal. 

■Attempts  to  determine  the  direction  and  velocity  of  LP  noise  in  Norway  by  frequency- 
wavenumber  methods  have  been  without  success  thus  far,  due  to  limited  data  coverage,  small 
number  of  sites  (three),  and  large  separations  (100  km). 

H.  W.  Briscoe 
R.  T.  Lacoss 


E.  LASA  STATION  CORRECTION  STUDIES 

We  have  not  reported  on  our  work  with  station  corrections  for  several  years,  although  they 
are  basic  to  the  operation  of  the  array  and  a  number  of  changes  have  been  made  in  the  way  they 
have  been  measured  and  used.  A  description  of  the  work  performed  in  this  reporting  period  will 
be  preceded  here  by  a  summary  of  previous  work  not  mentioned  in  preceding  Semiannual  Reports. 

When  the  arrival  times  of  P  at  each  LASA  subarray  are  measured,  for  a  well-recorded 
event,  and  compared  with  the  times  predicted  by  standard  tables  and  published  source  param¬ 
eters,  two  kinds  of  error  are  observed.  First,  the  arrival  time  at  any  given  subarray,  such 
as  the  "reference  subarray"  Aft,  differs  from  the  predicted  time,  and  second,  the  measured 
moveouts,  relative  to  this  reference,  differ  from  the  corresponding  computed  moveouts.  The 
first  error  contains  source  parameter  errors,  propagation  anomalies  in  the  vicinity  of  the  source 
(source  correction),  departure  of  actual  travel  time  from  world-wide  average  tables,  record 
reading  error,  and  the  absolute  station  correction,  relative  to  the  standard  tables  used,  of  the 
reference  subarray.  The  absolute  station  correction  is  presumably  due  to  local  propagation 
anomalies  in  the  vicinity  of  the  reference  site,  although  the  separation  of  "local"  effects  from 
"regional  travel  time  variations"  is  necessarily  imprecise.  The  second  error,  for  any  given 
subarray  other  than  A(^,  is  attributed  to  the  "relative  station  correction,"  although  all  the  sources 
of  error  enumerated  above  are  present  to  some  degree.  One  expects  that  source  parameter 
errors  and  source  corrections  are  comparatively  less  important  in  the  relative  station  correc¬ 
tion,  and  that  only  departures  of  travel  time  derivatives  from  standard  tables  will  be  relevant 
over  the  LASA  aperture.  However,  the  local  propagation  anomalies  under  the  subarray  in  ques¬ 
tion  and  under  the  reference  site  will  both  contribute. 

Earthquake  arrival  time  data,  yielding  information  on  both  kinds  of  station  correction,  is 
valuable  both  for  providing  insight  into  the  velocity  structure  under  the  array  and  for  the  opera¬ 
tional  determination  of  epicenters  by  either  triangulation  or  beamsplitting.  In  this  Section  we 
describe  our  work  with  this  arrival  time  data,  the  aim  of  which  is  to  reduce  it  to  a  form  more 
useful  for  both  purposes. 

Our  earliest  work  with  station  corrections  was  based  upon  measurements  of  P  arrival  times, 
taken  from  film  recordings  at  the  Laboratory,  and  reduced  to  relative  station  correction  data  by 
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the  use  of  USCGS  reports  on  source  parameters.  The  moveout  anomaly  data  was  displayed,  for 
each  subarray,  as  a  function  of  event  bearing  only  (since  the  data  depend  most  strongly  on  bearing) 
and  fit  by  least  squares  to  a  five-term  Fourier  series  in  the  bearing  angle.  Values  from  the 
resulting  curves  were  then  used  for  epicenter  determination  bv  triangulation.  It  was  evident  at 
once  that  the  results  had  iiiree  shortcomings;  (1)  the  data  base  of  about  ISO  events  was  inade¬ 
quate,  (2)  it  would  be  necessary  to  know  the  station  corrections  as  functions  of  both  distance  and 
bearing,  and  (3)  the  extremely  non-uniform  geographical  character  of  natural  seismicity  meant 
that  any  simple  curve-  or  surface-fitting  technique  would  lead  to  results  dominated  by  the  data 
from  a  small  fraction  of  the  teleseismic  epicentral  area. 

In  order  to  alleviate  the  first  two  shortcomings,  we  undertook  to  obtain  data  on  a  large 
population  of  clearly  recorded  events  and  to  compute  average  corrections  for  each  site  as  func¬ 
tions  of  distance  and  bearing.  This  effort  was  carried  out  under  subcontract  by  the  Seismic  Data 
Laboratory  (Teledyne)  and  resulted  in  time  and  amplituoe  anomaly  readings  for  over  600  well- 
recorded  events  (of  average  magnitude  5.1).  In  order  to  use  the  resulting  corrections  in  any 
location  programs,  it  was  necessary  either  to  fit  the  data  with  some  readily-computable  function 
of  two  variables,  o'"  store  the  numbers  on  some  fixed  grid  of  points  in  the  distance-bearing  plane 
and  use  a  table  lookup  procedure  when  corrections  are  needed.  For  the  automatic  event  location 
program  in  use  at  the  LASA  Data  Center,  it  was  decided  to  use  the  second  approach,  although 
PDP-7  memory  limitations  restricted  us  to  rather  coarse  cells.  We  used  5°  in  distance  and  30° 
in  bearing  to  define  216  distance-bearing  cells,  and  averaged  the  raw  anomaly  measurements 
over  each  cell  that  contained  data.  At  grid  points  corresponding  to  cells  containing  no  data,  the 
station  corrections  were  put  equal  to  zero. 

The  results  were  fairly  satisfactory  for  triangulation  on  events  which  fell  in  occupied 
cells,  but  errors  as  large  as  10°  or  more  could  occur  for  events  in  unoccupied  distance-bearing 
cells.  It  was  also  found  that  this  form  of  the  station  corrections  was  inadequate  when  patterns 
of  closely -packed  beams  were  formed  and  their  outputs  used  for  epicenter  determination,  as 
described  in  Sec.  II-A.  On  contour  plots  of  beam  output  energy  vs  latitude  and  longitude,  abrupt 
changes  at  cell  boundaries  can  be  seen,  especially  in  the  transition  from  occupied  to  unoccupied 
cells. 

Our  next  step  (at  the  beginning  of  the  present  reporting  period)  was  to  return  to  the  raw 
station  anomaly  data  (in  the  form  of  a  tape)  and  generate  histogram  displays  of  the  station  anom¬ 
aly  for  each  site.  Bounds  were  then  chosen  for  extreme-value  data  from  each  site  by  studying 
these  histograms,  and  the  data  arbitrarily  truncated  within  each  pair  of  bounds.  Next,  a  two- 
dimensional  sliding  window,  5°  by  5°  in  the  distance-bearing  plane  was  passed  over  the  data  to 
obtain  average  values  on  a  grid  of  1296  distance-bearing  pairs.  Grid  points  corresponding  to 
cells  containing  no  data  are  flagged  and  station  correction  numbers  for  these  cells  are  made  up 
in  the  following  way.  First,  a  value  is  assigned  to  each  empty  cell  equal  to  the  average  correc¬ 
tion  for  the  nearest  occupied  cell,  using  great -circle  distance  on  the  earth's  surface.  Next,  this 
complete  grid  of  values  is  passed  repeatedly  thr''ugh  a  two-dimensional  smoothing  program  which 
is  not  allowed  to  change  values  corresponding  to  occupied  cells  (i.e.,  real  data).  The  smoothing 
algorithm  is  Laplacian;  a  "new  value"  is  computed  for  each  grid  point  as  the  average  of  the  old 
value  and  that  of  its  four  nearest  values,  ar.d  the  smoothed  value  actually  uoud  is  a  weighted 
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Section  II 


aveiage  of  the  old  and  "new"  numbers.  Afte;;  many  iterations,  the  result  approaches  the  solution 
of  the  finite-mesh  I.aplace  equation  with  the  original  real  data  as  boundary  values.  A  moderate 
amount  of  smoothing  is  performed  at  the  end.  holding  no  points  fixed.  The  result  is  a  smooth 
lunction  oi  distance  and  bearing,  with  values  at  each  grid  point  and  still  very  close  to  the  original 
data  uhere  such  e.xisted.  A  contour  diagram  of  the  F4  correction,  using  distance  and  bearing  as 
polar  coordinate.^,  is  shown  in  h'ig.  17. 

The  new  coi-rections  have  so  far  been  ised  operationally  only  in  beamforming  programs, 
where  considerably  .smonthei-  patterns  are  now  found.  However,  these  interpolated  smoothed 
corrections  arc  also  useful  in  earth  structure  studies,  because  they  represent  the  data  in  a  way 
which  IS  not  strongly  biased  by  natural  seismicity.  One  such  application  is  the  display  of  move- 
out  diagrams.  Tbe.se  are  produced  by  choosing  an  epicenter  a  d  computing  station  corrections 
from  the  grid  valves  for  each  site  as  well  as  the  relative  distance  of  r  ch  site  to  the  epicenter. 
Alien  the  corrections  are  plotted  vs  projected  epicentral  distance,  the  result  is  a  "typical  move- 
out  piot,  which  is  representative  of  tlie  moveouis  obtained  for  events  near  that  epicenter  with 
the  theoretical  moveouts  removed.  The  example  sliown  in  Fig.  18,  for  an  epicenter  at  72.5*, 
bearing  H5'.  clearly  suggests  a  set  of  dipping  layers  for  the  L.ASA  substructure.  Attempts  to 
explain  ihe  station  correction  results  are  being  made,  based  on  models  of  the  structure  obtained 
from  such  movemit  plots. 

U  e  are  now  re-referoncing  station  correction  data  on  the  1968  Herrin-Taggart  travel  times 
and  applying  the  programs  described  above  to  the  absolute  station  correction  data  at  A0. 

B.J.  Kelly 


.j.  18.  "Moveoot  o.ogrom"  for  bearing  145®,  distance  72.5®. 
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m.  GEOPHYSICAL  DATA  HANDLING  TECHNIQUES 


A.  DATA  ANALYSIS  CONSOLE 

The  Data  Analysis  Console  described  in  the  previous  Semiannual  Technical  Summary  Report 
(31  December  1967)  has  had  two  programs  added  to  it  which  are  described  below.  These  addi¬ 
tions  complete  the  console  software  design  for  our  present  hardware  configuration.  A  complete 

14 

description  of  this  current  system  is  given  in  a  recent  report. 

Additional  hardware  (expanded  core  memory,  3-million-word  drum,  automatic  priority  in¬ 
terrupts  and  block  tran.,fer  scope)  has  been  ordered  to  greatly  increase  the  capability  of  the 
basic  PDP-7  computer.  A  new  software  package  will  be  designed  using  the  new  hardware  which 
will  create  a  new  data  analysis  console  that  will  be  similar  to  the  one  in  current  use,  but  will 
be  faster  and  allow  more  powerful  and  versatile  processing. 

The  two  new  features  of  the  present  console  are  as  follows. 


1.  Event  Detecting  {J] 

This  program  scans  the  data  tape  and  types  out  the  time  of  all  teleseisms  encountered. 
The  event  detection  algorithms  have  been  described  previously.^ ^ 


2.  Matched  Filtering  {M] 


This  program  cross  -orrelates  the  reference  trace  with  a  matched  filter  that  is  specially 
designed  to  enhance  the  SNR  of  long-period  dispersed  surface  waves.^  The  matched  filter  is  a 
sinusoid  with  a  linear  frequency  modulation  (chirp).  The  matched  filter  has  three  parai..eters 
that  can  be  input  to  the  program  via  the  telety-pewriter;  the  period  at  the  start  of  the  chirp,  the 
period  at  the  end  of  the  chirp,  and  the  length  of  vhe  chirp.  The  maximum  allowable  length  of  the 
chirp  matched  filter  impulse  response  is  830  samples  which  corresponds  to  160  to  1920  seconds 
depending  on  the  sampling  chosen  in  the  Initialize  program.  Enough  of  the  input  reference  will 
be  convolved  with  this  chirp  to  give  512  samples  of  filtered  output,  which  will  appear  in  the  ref¬ 


erence  buffer. 


P.  L.  Fleck 
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